#include "Mesh.h"
#include <iostream>
#include <Eigen/Dense>
#include "Quadrature.h"
#include "FEBase.h"
#include "InitialCondtion.h"
#include "Equation.h"

using Eigen::MatrixXd;
int main()
{

	Mesh mesh(10);
	Equation equation(mesh);
	equation.addVariable(SECOND);
	equation.addInitCondtion();
	equation.computeUlUr();
	std::cout << equation.postPoint() << endl;
	equation.outputExact("exact");
	equation.output("file");
	
	equation.solve(100);

}